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With the equations of state provided by the newly developed density dependent relativistic 
Hartrce-Fock (DDRHF) theory for hadronic matter, the properties of the static and (3- 
equilibrium neutron stars without hyperons are studied for the first time, and compared to 
the predictions of the relativistic mean field (RMF) models and recent observational data. 
The influences of Fock terms on properties of asymmetric nuclear matter at high densities 
are discussed in details. Because of the significant contributions from the a- and w-exchange 
terms to the symmetry energy, large proton fractions in neutron stars are predicted by the 
DDRHF calculations, which strongly affect the cooling process of the star. The critical 
mass about 1.45 AIq, close to the limit 1.5 Mq determined by the modern soft X-ray data 
analysis, is obtained by DDRHF with the effective interactions PK02 and PK03 for the 
occurrence of direct Urea process in neutron stars. The maximum masses of neutron stars 
given by the DDRHF calculations lie between 2.45 Mq and 2.49 Mq, which arc in reasonable 
agreement with high pulsar mass 2.08 ± 0.19 Mq from PSR B1516+02B. It is also found 
that the mass-radius relations of neutron stars determined by DDRHF arc consistent with 
the observational data from thermal radiation measurement in the isolated neutron star 
RX J1856, QPOs frequency limits in LMXBs 4U 0614-^09 and 4U 1636-536, and redshift 
determined in LMXBs EXO 0748-676. 

PACS numbers: 21.30.Fe, 21.60.Jz, 21.65.-f, 26.60.-c 



I. INTRODUCTION 



The investigative territory of nuclear physics has been enormously expanded with the con- 
struction of the new accelerator facihties as well as the development of the land- and space-based 
observatories. The exploration over the phase diagram of matter has been extended to the ex- 
treme conditions of density, pressure and temperature during the last several decades, which is 



2 



now one of the hottest topics in both theoretical and experimental nuclear physics. While around 
the saturation density, the nuclear matter properties can be well calibrated by terrestrial exper- 
iments with the atomic nuclei, the neutron stars are the natural laboratories in the unverse for 
exploring the equation of state (EoS) of baryonic matter at low temperature and higher baryonic 
densities [l, [3, [s]. In addition, the probe to the elliptical flow and kaon production in heavy- ion 
collisions provides extra information for the nuclear EoS at high temperature and about 2 ~ 4.5 



times nuclear saturation density 



aa. 



As one of the most exotic objects in the universe, neutron star plays the role of a bridge between 
nuclear physics and astrophysics. With the discovery of neutron in 1932, the concept of the neutron 
star were firstlyproposed by Landau Q]- Two years later, the neutron star was deemed to be formed 
in supernovae [s]. In 1960s, the observed radio pulsars were identified as the rotating neutron 
stars {lo| . Currently the neutron star is generally considered to be of the crust structure 11 1 . Below 
the atmosphere and envelope surface with a negligible amount of mass, the crust extends about 1 
to 2 km into star, which mainly consists of nuclei and free electrons. With the density increasing, 
the dominant nuclei in the crust vary from ^^Fe to extremely neutron-rich nuclei and neutrons 
may gradually leak out of nuclei to form the neutron fluid. The outer core {p > Po/3) of neutron 
stars is composed of a soup of nucleons, electrons, and muons. In the inner core, exotic particles 
may become abundant, such as the strangeness-bearing hyperons and/or Bose condensates (pions 
or kaons), and a transition to a mixed phase of hadronic and deconfined quark matter becomes 
possible. Although similar EoS at saturation and subsaturation densities is obtained by various 
nuclear matter models, their deviations are very remarkable in the high density region, which is 
very essential in describing and predicting the properties of neutron stars. Further investigations 
are therefore necessitated for the detailed structure over the density range of neutron stars. 

The recent observations of neutron stars have been reviewed, such as in Ref. [3]. The existence 
of massive compact stars of 2 Mq or above is now unveiled by some evidence. Careful analysis of the 
Rossi X-ray Timing Explorer {RXTE) data for the quasi-periodic brightness oscillations (QPOs) 
discovered from low-mass X-ray binaries (LMXBs) 4U 1636-536 shows that several neutron stars in 
LMXBs have the gravitational masses between 1.9 Mq and possibly 2.1 Mq [3]. Measurements on 
millisecond pulsars in globular cluster NGC 5904 (M5) during 19 years of Arecibo timing yield M = 
2.08 ± 0.19 Mq for PSR B1516-F02B [ij. Whereas a much larger pulsars mass 2.74 ± 0.21 Mq is 
presented very recently for PSR J1748-2021B in NGC 6440 [l^. Besides the maximum mass limits, 
the mass-radius relation is also constrained by the recent observations. The thermal radiation 
spectra in X-rays and in optical-UV from the isolated neutron star RX J1856. 5-3754 (shorthand: 
RX J1856) determine the large radiation radius Roo as 16.8 km [l^. The model fitting to the 
high-quality X-ray spectrum of the quiescent LMXB X7 in the globular cluster 47 Tuc prefers a 
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171]. In another LMXB, EXO 0748- 



rather large radius of 14.5l'lj^'g km for a 1.4 Mq compact star 
676, a pair of resonance scattering lines consistent with Fe XXV and XXVI, gives the redshift z 
about 0.345, which constrains the mass M > 2.10 =b 0.28 Mq and the radius R > 13.8 it 1.8 km for 
the same object [isl . In addition, the highest QPOs frequency 1330 Hz ever observed in 4U 
0614+09 implies the mass M < 1.8 Mq and the radius i? < 15 km in this object [2^. Furthermore, 



modern observational soft X-ray data of coohng neutron stars associated with popu 



models' analysis reveal that an acceptable EoS s 
in neutron stars with masses below 1.5 Mq 



na. 



22 



23 



ar synthesis 



2l| to occur 



1 not allow the direct Urea process 
[2^ . All of these indicate the strict constraints 



on the EoS of strongly interacting matter at high densities. 

For the description of nuclear matter and finite nuclei, the relativistic many-body theory has 
achieved great success during the past years. One of the most successful representatives is the 
relativistic Hartree approach with the no-sea approximation, namely the relativistic mean field 



(RMF) theory 



25 



26, 



271 ]. With a limited number of free parameters including the meson masses 



and meson-nucleon coupling constants, the appropriate quantitative de scr i ptio ns are o 
RMF for bot 
3: 
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33, 



1 stable nuclei and exotic ones with large neutron excess 
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4QI, 



ill. 



281, 



29 



3C 



31 



jtained by 
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33 



341, 



431 ] ■ a plenty of theo- 



After the first theoretical calculations of the neutron star properties |42 
retical prediction was made by both non-relativistic and relativistic approaches in the literature. 
At the early development of RMF, it was applied to evaluate the total mass and radius of neutron 



stars 



44| . In the further development, the nuclear medium effects were taken into account by 



introducing the explicit or implicit density dependence into the meson-nucleon couplings, i.e., the 



density dependent meson-nucleon couplings 



son fields 



39,Q,Q 



respectively. In Refs. [4 




471] and the non-linear self-couplings of the me- 
5l| , the effects of the non-linear self-coupling of 



a, uj and p mesons were studied in describing the nuclear matter and neutron stars. On the other 
side, the influence on mean field potentials, saturation properties of nuclear matter, the EoS, the 
maximum mass and radius of neutron stars was systematically investigated with explicit density 

53| . In addition, the consequences on compact sta r 



dependence in the meson-nucleon couplings 



52 



54 



55 



properties were studied with the inclusion of the degree of freedom of hyperons 
With more accurate experimental data of the neutron radius of ^'^^Pb, the correlation between the 
neutron skin thickness in finite nuclei and the symmetry energy of nuclear matter was discussed 



53, 



5S 



6OI ]. Further investigation in Ref. [61[ proposed the relation between the neutron skin 
thickness of a heavy nucleus and neutron star radius that the larger neutron skin thickness prefers 
the larger neutron star radius, which implies the constraints on the EoS and on the cooling mech- 
anism of neutron stars. Besides the RMF approach, the Dirac-Brueckner-Hartree-Fock (DBHF) 
and Brueckner-Hartree-Fock (BHF) with three-body force approaches were also applied to study 
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the neutron star properties with the reahstic nucleon-nucleon interactions 6^, IS^, ISJ, l65|, |66|] . 

In the RMF approach, however, the Fock terms are neglected by the reason of simphcity. 
From the recent development on the relativistic Hartree-Fock theory, i.e., the density dependent 
relativistic Hartree-Fock (DDRHF) theory 



67( 1 ■ it is found that the Fock terms are of special 
importance in determining the nuclear structure properties. Within DDRHF, the quantitatively 
comparable precision with RMF are obtained for the structure properties of nuclear matter and 
finite nuclei 



63, 



]. Particularly, the new constituents introduced with the Fock terms, i.e., 
the /9-tensor correlations and pion exchange potential, have brought significant improvement on 
the descriptions of the nuclear shell structures [68| and the evolutions [691]. Furthermore, the 



excitation properties and the non-energy weighted sum rules of the Gamow- Teller resonance and 
the spin-dipole resonance in the doubly magic nuclei have been well reproduced by RPA based on 
the DDRHF approach fully self-consistently [3]. Since the nuclear structure properties around the 
saturation density are evidently affected by the Fock terms, one might expect remarkable effects 
from the Fock terms on the nuclear matter properties in the high density region. Especially with the 
inclusion of the new ingredients in DDRHF, remarkable adjustment occurs on the coupling strength 
of the dominant mean fields (g^ and g^), which may bring significant effects when exploring to the 
high density region. 

In this paper, the properties of the static and /^-equilibrium neutron stars without hyperons 
are studied within the DDRHF theory. As compared to the calculations of the RMF theory, the 
applicable ranges of density and isospin asymmetry are tested for DDRHF as well as the consistence 
with recent observational constraints of compact stars. In Section |TI] is briefly introduced the 
formulism of DDRHF for nuclear matter and neutron star. In Section IIIH the calculated results 
and discussions are given, including the properties of symmetric and asymmetric nuclear matter 
in comparison with RMF in Section [HI Al in which the effects of Fock terms are studied in details, 
and the investigations of neutron stars in comparison with recent observational data in Section 
IIIIBI Finally a summary is given. 

II. GENERAL FORMULISM OF DDRHF IN NUCLEAR MATTER 



The relativistic Hartree-Fock (RHF) theory with density dependent meson-nucleon couplings. 



i.e., the DDRHF theory, was firstly introduced in Ref. 
effective interactions can be found in Refs. 



671, 



68 



6S 



6711. and the applications and corresponding 



7ll ]. In the following we just briefly recall 



the general formulism of DDRHF in nuclear matter and the application in neutron stars. More 
details of the RHF theory are referred to Refs. 



671, 



6g 



72| 



As the theoretical starting point, the Lagrangian density of DDRHF is constructed on the one- 
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boson exchange diagram of the NN interaction, which contains the degrees of freedom associated 
with the nucleon (ip), two isoscalar mesons (a and uj), two isovector mesons (vr and p), and the 



photon {A). Following the standard procedure in Ref. Ty], one can derive the Hamiltonian in 
nucleon space as, 

H= f d^xi>{-n-V + M)tl; + ^ I cfxdSJ2^Pix)ij{y)T4l,2)D^il,2)i;{y)i;{x), (1) 

where (p = a,uj, pjir and A, and denotes the propagators of mesons and photon. The interacting 
vertex in the Hamiltonian ([1]) reads as, 

r,(l,2) = -5,(1)5,(2), (2a) 

r.(l,2) = +g,(l)7,.(l)5.(2)7n2), (2b) 

r,(l,2) = +5p(l)7^(l)r(l) -5^(2)7^2)7^(2), (2c) 

r^(l, 2) = [/^f757 • V]i • [Af757 • , (2d) 

r^(l,2) = +jW-t,Mj>^{1-ts)]2. (2e) 

In current work, the /o-tensor correlations are not enclosed. In the above expressions and following 
context, the isovectors are denoted by arrows and the space vectors are in bold type. 

In general, the time component of the four-momentum carried by mesons are neglected on the 
level of the mean field approximation. This neglect has no consequence on the direct (Hartree) 
terms, while for the exchange (Fock) terms it amounts to neglect the retardation effects. The 
meson propagators are therefore of the Yukawa form, e.g., in the momentum representation, 

^</>(l,2)= (3) 

where the exchanging momentum q = P2 — Pi, and (j) = cr,uj, p and vr. 

For the description of nuclear matter, the coulomb field thus could be neglected, and the mo- 
mentum representation is generally adopted in the Hamiltonian. Due to time-reversal symmetry 
and rotational invariance, the self-energy S can be expressed as 

S(p) = ^sip) + 7oSo(p) + 7 • pSy(p), (4) 

where p is the unit vector along p, and the scalar component S5, time component Sq and space 
component of the vector potential are functions of the four-momentum p = {E{p) , p) of nucleon. 
With the general form of the self-energy, the Dirac equation in nuclear matter can be written as, 

(7 • p* + M*) uip, s, t) = ^oE*u{p, s, r), (5) 
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with the starred quantities, 



p* = p + pT.v{p), 

*2 



(6a) 
(6b) 
(6c) 



which obey the relativistic mass-energy relation E* = p* + M* . With this relationship, one can 



introduce the hatted quantities as 



- P 
P =— 



M =- 



M* 



E* 



(7) 



With the momentum representation, the Dirac equation ([5]) can be formally solved and the Dirac 
spinors with positive energy read as 

/ 



u{p,s,t) 



E* + M* 
2E* 



1/2 



1 

cr ■ p* 



XsXr, 



(8) 



\ E* + M* 

where Xs and Xt respectively denote the spin and isospin wave functions. The solution ([8]) is 
normalized as 



u'{p,s,t)u{p,s,t) = 1. 



(9) 



The stationary solutions of the Dirac equation ([5]) consist of the positive and negative energy 
ones, and one can expand the nucleon field operator -0 in terms of Dirac spinors. Within the mean 
field approximation, the contributions from the negative energy states are neglected, i.e., the no-sea 
approximation. The nucleon field operate ip is therefore expanded on the positive energy set as, 



P,S,T 

V'^(x)=^nt(p,.,r)e*^'-ct 



'P,S,T ■ 



(10a) 
(10b) 



p,S,T 



where Cp^s,T and Cp^s,T are the annihilation and creation operators. With the no-sea approximation. 



the trial Hartree-Fock ground state can be constructed as. 



(11) 



P,S,T 



where |0) is the vacuum state. The energy functional, i.e., the energy density in nuclear matter, is 
then obtained by taking the expectation of the Hamiltonian with respect to the ground state |$o) 
in a given volume O, 



e = i (<l>o| H |$o) ^ (T) + E (^^) = + E + ' 



(12) 
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where (p = a,(jj, p, vr, and 

ek = ^u{p,s,T){-f ■p + M)u{p,s,t), (13a) 

p,S,T 

^0 = ^ X] X] u{pi,Sl,Tl)u{p2,S2,T2)T^{l,2)-^u{p2,S2,T2)u{pi,Sl,Ti), (13b) 



Pl,si,ri P2,S2,T2 



= yZ yZ U{pi, Sl,Ti)u{p2, S2,T2)T^{1,2)^^ 7Tu{pi,Sl,Ti)u{p2,S2,T2). (13c) 

pi,Sl,Tip2,S2,T2 <P 

In the energy functional, £k denotes the kinetic energy density, and and respectively 
correspond to the direct (Hartree) and exchange (Fock) terms of the potential energy density. 
With the Dirac spinors in Eq. Q, one can obtain the contributions of the energy density from 
each channel. The kinetic energy density and the direct terms of the potential energy density 

can be written as, 

'■kF,i 



Bk = Yj ^ f ^'\'^dp{pP + MM^ , 



i=n,p 

1 ^2 

D _ , ^ 9uj 2 

2 mi 

Bp — +o — 2Pb'ii 

2 



(14a) 

(14b) 
(14c) 
(14d) 



where the scalar density ps, baryonic density pb and the third component pb3 read as. 



Ps 



Pb 



Pb3 



IT- 



i=n,p 



vr ./o 



p'dpM{p), 



E 



t.3 

3^2^ 



i=n,p 

1.3 1.3 

'^F,n _ '^F,p 

37r2 37r2 ' 



(15a) 
(15b) 
(15c) 



with the fermi momentum kp,! {i = n,p). 

TABLE 1: The terms A^, and in Eq. H 







b4p,p') 


C4p,p') 


a 




9'Mp,p') 


-2gl^,{p,p') 


UJ 




-Agieup.p') 


-4gl^^{p,p') 


P 




~4glep{p,p') 


-4gp*p(p,p') 






-f^e.{p,p') 


2Sr [(p2 +p'2) <i>4p,p')-pp'Q4p,p')] 

TT 



Compared to the simple form of direct terms of the potential energy density, the exchange terms 
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are much more complicated. In the isoscalar channels {(j) = cr,Lo), the expressions read as, 

£f = ^-r' I ppdpdp' [a^{p,p') + M{p)M{p')B4p,p') + P{p)P{p')C^{p,p')] . (16) 



T,r' 



For the isovector channels = p, vr), one just needs to replace the isospin factor drr' by (2 — drr') 
in the above expression. The details of the terms and in Eq. (jl6p are shown in Table [H 

where the functions @^{p,p') and ^(i){p,p') are defined as, 

1 Att^ rnl+{p + p'f 4vr2 

dndi}'^—^ = ^in ; % ^ ^e<^(p,pO, (17a) 
"^■0 + Q PP ^0 + (P ~ F) PP 

f ■b v Att"^ {p^ + p'^ + "1^ 1 47r2 

/ df^dJ]'^^ = ^ , , ^ 9^(p,pO - 2 U 2^ci,^(p,p'). (17b) 

J + PP I 2pp J 

^From the potential energy densities in Eqs. (|13bp and (jl3c|) . one can perform the following 
variation, 

and obtain the self-energy S(p) which includes the direct terms, 

= -^Ps, (19a) 
= +4p6 + 4p63, (19b) 



and the exchange terms. 



KsiP) = j M(pOp>'E{'^-'[^- + ^-](P,P') + (2-'^-')[^P + ^^](py)}, (20a) 

j P'dp' {Srr' [A. + AJ(py) + (2 - [A, + , (20b) 



KviP) = I P{p')p'dp' {^rr' [C. + + (2 - 5^r') \C p + • (20c) 

In DDRHF, the explicit density dependence is introduced into the meson-nucleon couplings, i.e., 
the coupling constants g^^ g^^, gp and /tt are functions of the baryonic density ph- In the isoscalar 
meson-nucleon coupling channels 

following form, 

9<p (pb) = g<p (po) U (x) , (21) 

where x = Pb/po, and po is the saturation density of nuclear matter, and the function reads as 

l + b^{x + d^f 

U y^) - "■4>— — } — TTTi • (^^) 
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In addition, five constraints, (1) = 1, (1) = (1), and (0) = 0, are introduced to reduce 
the number of free parameters. For the isovector channels, the exponential density dependence is 
adopted for gp and /^r 

9p{Pb) = 5p(0)e-«^^ (23) 
fniPb) = /.(0)e-«-^ (24) 

Due to the density dependence in meson-nucleon couplings, the additional contribution, i.e., 
the rearrangement term S/j, appears in the self-energy S. In nuclear matter, it can be written as, 

E ^^T.^I[^(pKsip) + Koip)+Hp)Kvip)]p'dp. (25) 

^From Eqs. (119^ . (I20|) and (I25p , the scalar component S^, time component Eg and space 
component of the vector potential in Eq. ^ can be obtained as, 

^s{p) = + Sf(p), (26a) 
So(p) = S^ + Sf(p) + SH, (26b) 
Sy(p) = S^(p), (26c) 

from which the starred quantities in Eq. ([6]) and the hatted quantities in Eq. ([7]) can be obtained. 
Therefore, for nuclear matter with given baryonic density pb and neutron-proton ratio N/Z, one 
can proceed the self-consistent iteration to investigate their properties: with the trial self-energies, 
one can determine the starred quantities, hatted quantities and calculate the scalar density, and 
then get the new self-energies for next iteration until the final convergence. 

In this work, a neutron star is described as the /3-stable nuclear matter system, which consists 
of not only neutrons and protons, but also leptons A (mainly e~ and /x^). The equations of motion 
for the leptons are the free Dirac equations and their densities can be expressed in terms of their 
corresponding Fermi momenta, px = kp ■^/{3tt'^) (A = e~,/i~). The chemical potentials of nucleons 
and leptons satisfy the equilibrium conditions 

P-p =Pn - Pe, Pp =Pe, (27) 

where the chemical potentials pn, Pp, Pp and pe are determined by the relativistic energy- 
momentum relation at the momentum p = kp, 

Pi=^Q{kF,i) + E*{kF,i), (28a) 



where i = n,p and A = e~ , p~ . The lepton masses are respectively, rrie = 0.511 MeV and = 
105.658 MeV. In addition, the baryon density conservation and charge neutrality are imposed as 

Pb =Pn + Pp, Pp =Pp + Pe- (29) 
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With these constraints, the energy density of neutron stars is then obtained as 



(30) 



Here the leptons are treated as the free Fermi gas by assuming that there are no interactions 
between leptons and nucleons or mesons and the kinetic energies of leptons can be expressed as, 



1 

^ -10 



(31) 



With the thermodynamic relation, the pressure of the neutron star system can be obtained as 

d ens 



P{Pb) = Pt 



dpb Pb 



E 



Pi Pi - £n 



(32) 



n 



At low density region {pb < 0.08 fm-3), instead of DDRHF calculations, the BPS[7| and BBP[7| 
models are chosen to provide the proper EoS. 

The structure equations of a static, spherically symmetric, relativistic star are the Tolman- 



Oppenheimer-Volkov (TOV) equations 



42 



43l |. Taking c = G = 1, the TOV equations read as 



dP _ [P{r) + e{r)] [M(r) + 4Trr^P{r)] 
dr 
dM 



r[r- 2M (r)] 



dr 



-Attt e{r), 



(33a) 
(33b) 



where P{r) is the pressure of the star at radius r, and M(r) is the total star mass inside the sphere 
of radius r. Taking the equation of state of stellar matter as the input, one could proceed with the 
solution of TOV equations. The point R, at which the pressure vanishes, i.e., P{R) = 0, defines 
the radius of the star and the corresponding M{R) is the gravitational mass. For a given EoS, 
the TOV equation has the unique solution which depends on a single parameter characterizing the 
conditions of matter at the center, such as the central density p(0) or the central pressure -P(O). 

III. RESULTS AND DISCUSSIONS 



In this paper, the EoS and the neutron star properties are studied in DDRHF with the effective 
interactions PKOl, PK02 and PK03 



67 



69( 1 . As shown in Table HH the coupling constant 5^p(0) 



is fixed to the value in the free space in PKOl whereas free to be adjusted in PK02 and PK03, 
and vr-coupling is not included in PK02. For comparison, the results calculated by RMF are also 



discussed. The effective interactions used in RMF ca^ 
ones GL-97 144|, NLl [75^, NL3 [3^, NLSH |Y6|, TMl 



ones TW99 [38|, DD-MEl [40], DD-ME2 [4ll] and PKDD [39(]. 



culations include the nonlinear self-coupling 



4g|] anc 



PKl 



39|, and the density-dependent 
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TABLE II: The effective interactions PKOl, PK02 and PK03 of DDRHF[67|l, where the masses M = 938.9 
MeV, = 783.0 MeV, rup = 769.0 MeV and = 138.0 MeV. 







9a 


9ui 


<?p(0) 


/x(0) 


flp 




Po 


PKOl 
PK02 
PK03 


525.769084 
534.461766 
525.667686 


8.833239 
8.920597 
8.895635 


10.729933 
10.550553 
10.802690 


2.629000 
4.068299 
3.832480 


1.000000 
1.000000 


0.076760 
0.631605 
0.635336 


1.231976 
0.934122 


0.151989 
0.151021 
0.153006 




la 


ha 




da 








du, 


PKOl 
PK02 
PK03 


1.384494 
1.375772 
1.244635 


1.513190 
2.064391 
1.566659 


2.296615 
3.052417 
2.074581 


0.380974 
0.330459 
0.400843 


1.403347 
1.451420 
1.245714 


2.008719 
3.574373 
1.645754 


3.046686 
5.478373 
2.177077 


0.330770 
0.246668 
0.391293 



A. Properties of nuclear matter 

1. Bulk properties 

In Table HIH are shown the bulk quantities of nuclear matter at saturation point, i.e., the sat- 
uration density the binding energy per particle Eb/A, the incompressibility K, the symmetry 
energy J and the scalar mass Mg/M. The results calculated by RMF with both the nonlinear 
self-coupling effective interactions and the density-dependent ones, which have been studied sys- 
tematically in Ref. 



53( 1 ■ are included for comparison. The saturation density and the binding 
energy per particle given by DDRHF with PKO series are around 0.152 fm"^ ^nd -16.0 MeV, 
respectively, close to the values provided by RMF. The incompressibility K calculated by DDRHF 
with PKOl, PK02 and PK03 range from 249 to 262 MeV, close to the values given by RMF with 
density dependent effective interaction. In contrast, relatively large values of K (270~360 MeV) 
are obtained by RMF with the non-linear self-coupling of mesons, except GL-97 and NLl. For 
the symmetry energy J, the non-linear version of RMF also presents relatively large values (36~44 
MeV) except GL-97, whereas the density-dependent version of RMF (except PKDD) provides com- 
parative values (32~34 MeV) to DDRHF with PKOl, PK02 and PK03. For the scalar mass M^, 
GL-97 gives the largest value and TW99 presents the smallest. The values given by DDRHF with 
PKO series are around 0.60, close to those by RMF with the non-linear self-couplings of mesons 
except GL-97, and systematically smaller values are obtained by RMF with the density dependent 
meson- nucleon couplings. 
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TABLE III: The saturation density po (fni^"^), binding energy per particle Eb/A (McV), incompressibility 
K (McV), asymmetry energy coefficient J (MeV) and the scalar mass Mg/M of nuclear matter at saturation 
point. 





PO 


Eb/A 


K 


J 




PKOl 


0.1520 


-15.996 


250.239 


34.371 


0.5900 


PK02 


0.1510 


-16.027 


249.597 


32.492 


0.6025 


PK03 


0.1530 


-16.041 


262.469 


32.987 


0.5862 


GL-97 


0.1531 


-16.316 


240.050 


32.500 


0.7802 


NLl 


0.1518 


-16.426 


211.153 


43.467 


0.5728 


NL3 


0.1483 


-16.249 


271.730 


37.416 


0.5950 


NLSH 


0.1459 


-16.328 


354.924 


36.100 


0.5973 


TMl 


0.1452 


-16.263 


281.162 


36.892 


0.6344 


PKl 


0.1482 


-16.268 


282.694 


37.642 


0.6055 


TW99 


0.1530 


-16.247 


240.276 


32.767 


0.5549 


DD-MEl 


0.1520 


-16.201 


244.719 


33.065 


0.5780 


DD-ME2 


0.1518 


-16.105 


250.296 


32.271 


0.5722 


PKDD 


0.1496 


-16.268 


262.192 


36.790 


0.5712 



2. Density dependence of the coupling constants 

In DDRHF, the medium effects are evaluated by the density dependence in the meson-nucleon 
couplings. To understand the EoS, it is worthwhile to have a look at the density dependence of 
the coupling constants. In Fig. [T]are shown the coupling constants ga, g^^, gp and fj^ as functions 
of baryonic density pb, where the results of the DDRHF effective interactions PKOl, PK02 and 
PK03 are given as compared to the RMF ones TW99, DD-ME2 and PKDD. As seen from Fig. [U 
all the effective interactions present strong density dependence in the low density region (p^ < 0.2 
fm~^) for both isoscalar (a and uj) and isovector (p and vr) meson-nucleon couplings. When density 
grows higher, g^ and g^^ in the left panels become stable. While due to the exponential density 
dependence, the isovector ones gp and fj^ tend to vanish (except gp in PKDD and PKOl) as shown 
in the right panels. From this aspect, one can understand that the isoscalar mesons provide the 
dominant contributions in the high density region. Compared to the RMF effective interactions, 
PKOl, PK02 and PK03 have smaller g^j and g^. This is mainly due to the effects of Fock terms, 
which lead to the recombination of the ingredients in the nuclear interactions. With the inclusion 
of Fock terms, the nuclear attractions are shared by the Hartree terms of cr-coupling and the 
Fock terms of to-, p- and 7r-couplings, and the repulsions are contributed by the Hartree terms of 
(j-coupling and the Fock terms of cr-coupling. While in RMF, the attraction and repulsion are 
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provided only by the Hartree terms of the a- and w-couphngs, respectively. 




0.0 0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.0 0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 

Pb (fm"') (fm"') 

FIG. 1: (Color online) The coupling constants g^, g^^, gp and /^r as functions of the baryonic density pf, 
(fm-3) for the DDRHF effective interactions PKOl, PK02 and PK03, and RMF ones PKDD, TW99 and 
DD-ME2. The shadowed area represents the empirical saturation region pb ~ 0.166 ± 0.018 fm^'''. 

It is not enough to adjust the isospin properties only within the nuclear saturation region. It 
is expected that the investigations on the EoS at higher densities and neutron star properties 
could provide the additional constraint. For the isovector coupling constants in the right panels 
of Fig. [H PKOl and PKDD present slightly weak density dependence in gp because of the fairly 
small density dependent parameter Up. In analogy to g„ and gi_j, the RMF effective interactions 
give larger values of gp. This is also due to the exchange contributions. In DDRHF, significant 
contributions to the isospin part of nuclear interaction are found in exchange terms of isovector 
mesons as well as isoscalar ones. It is different from the situation in RMF that the isospin properties 
are only described by the direct part of p-coupling. For the vr-meson, the contribution in neutron 
stars is negligible since fn tends to vanish at high densities. 



3. Equations of state 

The equations of state calculated by DDRHF with PKOl, PK02 and PK03 are shown in 
Fig. [2] and Fig. [3l respectively for the symmetric nuclear matter and pure neutron matter. The 
results calculated by RMF with TW99, DD-ME2 and PKDD are also shown for comparison. It 
is recommended to see Ref. [ssl for the density dependence of the EoS on more RMF effective 
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interactions. Seen from these two figures, identical behaviors of the EoS are provided by all the 
effective interactions at low density region [p^, < po) while in high density region exist pronounced 
deviations among different effective interactions. 
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FIG. 2: (Color online) The binding energy per particle Eb/A as function of the baryonic density pt, for sym- 
metric nuclear matter. The results are calculated by DDRHF with PKOl, PK02 and PK03, in comparison 
with those by RMF with TW99, DD-ME2 and PKDD. 



For the symmetric nuclear matter in Fig. [21 DDRHF with PKOl, PK02 and PK03 provides 
similar EoS as RMF with PKDD and DD-ME2, whereas much softer EoS is obtained by RMF 
with TW99 when density grows high. For the pure neutron matter in Fig. [3l the curves can be 
classified into three groups according to the behaviors of the EoS at high density region. Among 
all the effective interactions, the DDRHF ones present the hardest equations of state and the RMF 
one TW99 gives the softest one, whereas DD-ME2 and PKDD provide similar equations of states, 
which lie between the hardest and softest. Since the DDRHF parameterizations were performed 
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by fitting the properties of finite nuclei and nuclear matter around saturation point 
corresponds to the low density region, it becomes necessary to test the extrapolation of the effective 
interactions PKOl, PK02 and PK03 to high densities. 



4-. Symmetry energy 

The EoS property of isospin asymmetric nuclear matter is still ambiguous more or less. Different 
theoretical models predict quite different behaviors of the EoS for pure neutron matter. In most 
cases, it is due to the effective interactions obtained by fitting the properties of doubly magic 
nuclei, which have an isospin close to that of the symmetric nuclear matter. From this point of 
view, it becomes necessary to introduce the constraints, either from isospin asymmetric heavy-ion 
collisions experiments or from the data of nuclei with extreme isospin, into the fitting procedures 
of the effective interactions. 
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FIG. 3: (Color online) Similar to Fig. [2] but for pure neutron matter. 

The symmetry energy is an important quantity to illustrate the property of asymmetric nuclear 
matter. In general, the energy per particle of asymmetric nuclear matter E{pf,, (3) can be expanded 
in a Taylor series with respect to /?, 

E{p^, (3) = Eoipk) + (i'Esipb) + • • • , (34) 

where (3=1 — Ip^j p\, is the asymmetry parameter depending on the proton fraction. The function 
E(){ph) is the binding energy per particle in symmetric nuclear matter, and the symmetry energy 
Es{Pb) {J = Es{po)) is denoted as 

ld^E{pb,/3) 



Es{pb 



(35) 



/3=0 



2 5/32 

The empirical parabolic law in Eq. (j34|) is confirmed to be reasonable in all the range of the 
asymmetry parameter values, while at high density deviation from such a behavior is found [77 1. 



Fig. [3] shows the symmetry energy as a function of the baryon density p;,. The results are calcu- 
lated by DDRHF with PKOl, PK02 and PK03, in comparison with those by RMF with TW99, 
DD-ME2 and PKDD. As shown in Fig. d both DDRHF and RMF effective interactions present 
identical behaviors of the symmetry energy at low densities {p < pq), while sizeable enhancements 
in high density region are obtained by DDRHF with PKOl, PK02 and PK03 as compared to the 
RMF results. Among the RMF calculations, PKDD shows harder behavior than DD-ME2 and 
TW99, which provide identical symmetry energy in the whole density region. 

^From the energy functional in nuclear matter in Eq. (jl2p . one can obtain the contributions 
from different channels to the symmetry energy Es as. 



Es = Es,k + ^ {Es^^ + E^ 



E \ 



(36) 



where (p = a,UJ, p and vr. In fact, the direct terms of u-meson coupling have no contribution to the 
symmetry energy because of the nature of isoscalar-vector coupling. It is also expected that the 
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FIG. 4: (Color online) The nuclear symmetry energy Es (MeV) as a function of the baryon density pt, 
(fm"3). The results are calculated by DDRHF with PKOl, PK02 and PK03, in comparison to those by 
RMF with TW99, DD-ME2 and PKDD. 

one-pion exchange has minor effects since nuclear matter is a spin-saturated system. In Fig. [Sj the 
contributions from different channels to the symmetry energy are shown as functions of the baryon 
density pb- In the left panel are presented the contributions from the kinetic part and isoscalar 
channels, and only the results calculated by DDRHF with PKOl are shown in comparison with 
those by RMF with PKDD and DD-ME2. The contributions from the />-meson coupling are shown 
in the right panel, including the results calculated by DDRHF with PKOl, PK02 and PK03, and 
RMF with PKDD, DD-ME2 and TW99. 
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FIG. 5: (Color online) Contributions from different channels to the symmetry energy as a function of the 
baryon density pb- Left panel gives the contributions from the kinetic energy and isoscalar channels, and 
the ones from p-meson are shown in the right panel. See text for details. 

Within RMF, one can find that the kinetic part and the direct terms of a coupling, i.e., Es^k 
and Eg^, provide the dominant contributions to the symmetry energy. As shown in the left panel 
of Fig. [5l PKDD and DD-ME2 give identical values of Eg^^ + -^So-' ^^^'^ deviation between 
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them appearing in Fig. [J] is mainly due to their contributions from the p-couphng as seen from the 
right panel of Fig. [5l In the left panel of Fig. [H the values of Es^k + Eg^ given by the DDRHF 
calculations are found smaller than those of RMF with PKDD and DD-ME2. It is also seen that 
the Fock terms of a- and tij-couplings present significant contributions to the symmetry energy, 
which well interprets the stronger density dependence predicted by DDRHF than by RMF at the 
high density region (see Fig. [4]). As seen from the left panel of Fig. [5l the values of Eg^ increase 
rapidly in the low density region and tend to be stable about 20 ~ 25MeV at high density. While 
the exchange terms of L<j-coupling provide negative contributions to the symmetry energy at low 
density and reach the minimum about —11 MeV at = 0.1 fm~'^. After that the values of E^^ 
grow up and become comparable to the values of Eg^^ + ^Sa several times of the saturation 
density. 

In the right panel of Fig. [H the contributions of p-coupling, i.e., Es^p = Eg^ + Eg are found 
to be important for the symmetry energy at low density region. When density goes high, the values 
of Es^p given by all the effective interactions except PKOl and PKDD tend to zero due to their 
strong exponential density dependence in /j-nucleon coupling. Because of much smaller value of ap, 
PKOl presents larger contributions than PK02 and PK03 and contributes a value about 10 ~ 15 
MeV in the high density region. Due to the same reason, PKDD also provides larger values of 
Es^p than other two RMF effective interactions, and it reaches the maximum about 26 MeV at 
Pb — 0.41 fm~^, then falls down slowly. Comparing the values of Es^p given by PK02 and PK03 
to those by DD-ME2 and TW99, the contribution of p-coupling Es^p is depressed systematically 
in DDRHF. Such kind of depressions also exist between the results of PKOl and PKDD, which 
are of similar density dependence in p-nucleon coupling. This could be understood from the fact 
that smaller values of gp are obtained with the inclusion of Fock terms as seen in Fig. [TJ In fact, 
not only the p-meson but all the mesons take part to the isospin properties and are in charge of 
producing the symmetry energy via the Fock channel. 

Concluding the above discussions, one can find that the Fock terms play an important role in 
determining the density dependent behavior of the symmetry energy. It is then expected that the 
important constraints on the symmetry energy and the EoS of asymmetric nuclear matter could 
be obtained from the study of neutron stars. 

B. Properties of neutron star 

In this work, the static and /^-equilibrium assumptions are imposed for the description of neutron 
stars. With the density increasing, the high momentum neutrons will /? decay into protons and 
electrons, i.e., ?i <-> p + e~ + 17^, until the chemical potentials satisfy the equilibrium p,p = pn — fJ-e- 
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When the chemical potential of electron //e reaches the limit of the muon mass, the lepton /i~ will 
appear. The reaction e~ ^ + + i^g implies the equilibrium between the e~ and fi~ chemical 
potentials, i.e., /ie = /^/i- 



1. Density distribution 

To keep the equilibrium among the particle chemical potentials, protons, electrons and muons 
will appear with the density increasing in neutron stars. In Fig. [6] are shown the neutron, proton, 
electron and muon densities in neutron stars as functions of the baryon density. The results 
are calculated by DDRHF with PKOl, PK02 and PK03, in comparison to those by RMF with 
TW99, DD-ME2 and PKDD. The density distribution of various components in RMF with both 
the nonlinear self-coupling effective interactions and the density-dependent ones have been studied 



systematically in Ref. 53]- As seen from Fig. [6l the thresholds of /x^ occurrence predicted by 
different effective interactions are very close to each another, roughly around ph = 0.12 fm~^. It 
is shown that all the densities keep increasing monotonously with respect to the baryonic density 
Pb- Similar as the situation in the EoS of nuclear matter, different effective interactions present 
identical trends at low densities (p^ < Po), while remarkable deviations exist in high density region 
between the DDRHF and RMF predictions. Seen from Fig. [6l the results given by different effective 
interactions can be classified into three groups, the DDRHF ones, PKDD, and the RMF ones DD- 
ME2 and TW99. For the proton, electron and muon densities, the strongest density dependence is 
presented by DDRHF with PKOl, PK02 and PK03, while the softest behaviors are provided by 
RMF with DD-ME2 and TW99. In contrast, the softest behavior on neutron density is predicted 
by the DDRHF effective interactions, whereas TW99 and DD-ME2 present the hardest. This kind 
of reversion can be well understood from the relations in Eq. (|29p among the densities. 

It is known that the density fractions of each components in neutron stars are rather sensitive to 
the symmetry energy, as illustrated by associating Fig.[6]with Fig. [4l Due to the strong effects from 
the exchange terms of w-coupling in high density region (see Fig. [5]), DDRHF with PKOl, PK02 
and PK03 shows stronger density dependence on the symmetry energy and then are obtained the 
harder behaviors on the proton, electron and muon density distributions, as compared to the RMF 
calculations. For the deviations between different effective interactions within one theoretical 
model, e.g., between PKDD and DD-ME2, they are mainly due to the p-coupling as shown in 
the right panel of Fig. [Sj where the p-meson coupling of PKDD shows larger contributions to 
the symmetry energy. As a conclusion, the harder behavior on the symmetry energy at high 
densities, more difficult the system becomes asymmetric and more easier neutrons decay into 
protons and electrons, which leads to smaller neutron abundance and larger proton, electron and 
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FIG. 6: (Color online) The neutron (up-left panel), proton (up-right panel), electron (lower-left panel) and 
niuon (lower-right panel) densities in neiftron star matter as functions of the baryon density p^, (fni~^). The 
results are calculated by DDRHF with PKOl, PK02 and PK03, in comparison to those by RMF with 
TW99, DD-ME2 and PKDD. 

muon abundances in neutron stars. 



2. Proton fraction and direct Urea constraint 

^From the density distributions in Fig. [6l one can extract the proton fraction x = pp/ {pp + pn) 
within the range of density of neutron stars. Fig. [7] shows the proton fraction x as a function 
of baryonic density p^, where the results calculated by DDRHF with PKO series are presented in 
comparison to those by RMF with TW99, DD-ME2 and PKDD. Due to the stiff behavior on the 
symmetry energy (see Fig. U]), stronger density dependence of the proton fraction x in neutron 
star matter is obtained by DDRHF than RMF as shown in Fig. [71 

The cooling mechanism of neutron stars, which is sensitive to the proton fraction, could bring 
significant information of asymmetric nuclear EoS. Direct Urea (DU) processes n — > p+e~ + i>e and 
p+e~ n+Vf. lead the star to cool off rapidly by emitting the thermal neutrinos. The threshold of 

DU 



y found as 11.1% ^ x^^ ^ 14 



the proton fraction x^^ for the DU process occurring can be easi 

with the momentum conservation and charge neutrality {21! . [24]. Seen From Fig. [71 the critical 
density p^^ for the DU process occurring depends on the EoS. Once the critical density p^^ is 
reached in the center of a neutron star for a given EoS, the star will be efficiently cooled via the DU 
process. It is found that the values of x^^ given by DDRHF calculations correspond to fairly low 
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FIG. 7: (Color online) Proton fractions x = Pp/{pp + Pn) in neutron star matter for different DDRHF and 
DDRMF effective interactions. The dotted line labeled with x^^ is threshold for happening direct Urea 
process. Here x°^ = 14.8% is taken by assuming muons in the massless limit. 



critical densities while the results calculated by RMF with TW99 and DD-ME2 do not support the 
DU process occurring at all. The DU critical star masses M^^ and corresponding central densities 
p^^(O) are marked in Fig. [9] by filled squares. 



According to the analysis in Refs. 
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231, 



2J] , if the DU process is taken as a possible mechanism 
for neutron star cooling, an acceptable EoS shall not allow it to occur in neutron stars with masses 
below 1.5 Mq, otherwise it will be in disagreement with modern observational soft X-ray data 
in the temperature-age diagram. As a weaker constraint, the limit > 1.35 Mq could be 

applied. From the mass limit M^^ , are then obtained the constraint over the EoS that the density 
dependence of the symmetry energy should not be too strong, and probably not too weak, either. 
In Table ITVl are given the critical neutron star mass M^^ and central densities p^^(O) from the 
DDRHF and RMF calculations, which support the occurrence of the DU cooling process in stars. 

Seen from Table ITVl rather small mass limits M^^ are obtained by RMF with the non-linear 
self-coupling of mesons while the DDRHF calculations with PK02 and PK03 provide larger values 
of M^^, which are very close to the limit 1.5 Mq mentioned above and satisfy the weak constraint 
that > 1.35 Mq. For the calculation with PKOl, the DU cooling process will occur at the 

fairly low mass 1.20 Mq and central density p^^^ ~ 0.28 fm""^, which can be interpreted by the 
contributions of the /j-meson coupling to the symmetry energy. For the Es^p in the right panel of 
Fig. [5l the p- meson coupling in PKOl still has remarkable effects in the high density region due to 
the weak density dependence of Qp (see Fig. [T|). Due to the same reason, the RMF calculation with 
PKDD also support the DU cooling process to occur at a low mass limit 1.26 Mq. In contrast, 
as seen in Fig. Q the occurrence of the DU cooling process is not supported at all by the RMF 
calculations with TW99 and DD-ME2 as well as DD-MEl. It is expected that the occurrence of 
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the DU process could be introduced as a possible constraint in the future parameterizations of 
both DDRHF and RMF, e.g., for the /9-meson coupling ((?p(0) and Op). 

TABLE IV: Critical neutron star masses and central densities p^^ (0) for the occurrence of the DU 

cooling process and the criterion of the DU constraint given by both DDRHF and RMF elfectivc interactions. 
Fulhlhnent (violation) of a constraint is indicated with +(— ). 





PKOl 


PK02 


PK03 


GL-97 


NLl 


NL3 


NLSH 


TMl 


PKl 


PKDD 




1.20 


1.45 


1.43 


1.10 


0.75 


1.01 


1.20 


0.96 


0.94 


1.26 


pDU(o) [fni-3] 


0.28 


0.33 


0.31 


0.33 


0.20 


0.23 


0.24 


0.24 


0.23 


0.33 


^ 1.5 Mq 
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3. Pressure and maximum mass of neutron star 



In Fig. [HI the pressures of neutron star matter calculated by DDRHF effective interactions are 
shown as functions of the baryonic density p^. The results with RMF ones GL-97, NL3, TW99 
in Ref. 



53( 1 and also PKl, DD-ME2, PKDD have been included for comparison. It is found that 
PKOl, PK02 and PK03 provide idential behaviors with each another over the density dependence 
of the pressure, which are also close to the behaviors predicted by RMF with PKDD and DD-ME2. 
Among all the DDRHF and RMF calculations, NL3 provides the strongest density dependence and 
the softest are presented by GL-97. The behaviors given by PKl and TW99 lie between the results 
of DDRHF with PKO series and RMF with GL-97. 

800 




0.6 

P. (fm ') 



FIG. 8: (Color online) The pressure of neutron star matter as a function of the baryon density pb (fm '^). 
The resuhs are calculated by DDRHF with PKOl, PK02, and PK03, in comparison to those by RMF with 
GL-97, NL3, PKl, TW99, DD-ME2, and PKDD. 



The variation of the pressure with respect to density is essential to understand the structure of 
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neutron stars. Stronger density dependence of the pressure at high densities would lead to larger 
value of maximum mass for neutron stars that can be sustained against collapse. In Fig. [U the 
neutron star masses calculated by DDRHF with PKOl, PK02, and PK03 are shown as functions 
of the central density p(0). For comparison, are also shown the results calculated by RMF with 



53| for more studies 



GL-97, NL3, PKl, TW99, DD-ME2, and PKDD, and one could refer to Ref. 
with a variety of RMF effective interactions. From Fig. [U it is found that the maximum masses 
given by the DDRHF calculations lie between 2.4 Mq and 2.5 Mq with the central densities around 
0.80 fm~^, which are close to the prediction of RMF with DD-ME2. Notice that these values are 
also compatible to the observational constraint (M = 2.08 ±0.19 Mq) from PSR B1516+02B Q. 
In Table |V] are shown the maximum mass limits M^ax and the corresponding central densities 
Pmax(O) extracted from Fig. [9l As consistent with the description of the pressure, the non- linear 
RMF effective interaction NL3 presents a rather large value of the maximum mass Mmax = 2.78 Mq 
with small central density yOmax(O) = 0.67 fm~^, while the smallest Mmax and the largest /9max(0) 
are obtained by RMF with GL-97 and TW99, which gives the softest behaviors of the pressure 
(see Fig. [8]). Seen from Table IVl the values of M^ax given by all the effective interactions are in 
the appropriate agreements with the constraint on the maximum mass from PSR B1516+02B. 
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FIG. 9: (Color online) Neutron star mass as a function of the central density for different DDRHF and RMF 
effective interactions. Filled stars denote the maximum mass configurations, filled squares mark the critical 
mass and central density values p^^ (Q) where the DU cooling process becomes possible. The light 

grey horizontal bands around 2.08 Mq denote the Icr confidence level for the mass measurement of PSR 
B1516+02B [l^. The mass region of typical neutron stars is between 1.0 A4q and 1.5 AIq. 



4-. Mass-Radius relation and observational constraint 

Recent astronomic observations also provide the constraints on the mass-radius relation of neu- 
tron stars. In this paper, four typical observations are adopted to test the theoretical calculations. 
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TABLE V: Maximum mass limits A'/max {Mq), the corresponding central densities Pmax(O) (fm^"^) and radii 
i?(Mniax) (km) for neutron stars calculated by DDRHF and RMF effective interactions. The radii (km) for 
1.4 Mq neutron stars arc shown as well. 





PKOl PK02 PK03 


GL-97 NLl NL3 NLSH TMl PKl 


TW99 DD-MEl DD-ME2 PKDD 


-^Anax 
Pmax(O) 
i?(M,„ax) 


2.45 2.45 2.49 
0.80 0.81 0.78 
12.4 12.3 12.5 


2.02 2.81 2.78 2.80 2.18 2.32 
1.09 0.66 0.67 0.65 0.85 0.80 
10.9 13.4 13.3 13.5 12.4 12.7 


2.08 2.45 2.49 2.33 
1.10 0.84 0.82 0.89 
10.7 11.9 12.1 11.8 


R{1AMq) 


14.1 13.8 13.9 


13.3 14.7 14.7 14.9 14.4 14.5 


12.4 13.2 13.3 13.7 



1. The large radiation radius = 16.8 km [Roo = Rj — 2GM / Rc^) from the isolated 
neutron star RX J1856 Q. 

2. The redshift z ~ 0.345, the mass M > 2.10 ± 0.28 M© and the radius R > 13.8 ± 1.8 km 
constraints in LMXBs EXO 0748-676 

3. M < 1.8 Mq and i? < 15 km constraints from the highest frequency of QPOs 1330 Hz ever 
observed in 4U 0614+09 Q- 

4. Several neutron stars in LMXBs have gravitational masses between 1.9 Mq and possibly 
2.1 Mq from the QPOs data analysis in LMXBs 4U 1636-536 Q. 

In Fig. [10] are shown the mass-radius relations of neutron stars calculated by DDRHF with 
PKOl, PK02 and PK03, and RMF with GL-97, NL3, PKl, TW99, DD-ME2 and PKDD. The 
results with more RMF effective interactions have been investigated in Ref. 
the selected observational constraints are marked with different colors and grids as shwon in Fig. 
[TOl The causality hmit that y^dp/de < 1 results in R > 2.9GM/c^ a , [3] and the corresponding 
region in Fig. [10] is maked in black. Compared to all the observational limits, it is found that 
better agreements are obtained by the DDRHF effective interactions than the RMF ones. Among 
the RMF results, GL-97 is excluded by the limits from RX J 1856, and TW99 is excluded by the 
limits from both RX J1856 and EXO 0748-676, while NL3 could not fulfil the constraint from 4U 
0614+09. If upper mass limit 2.1 Mq is taken in 4U 1636-536, GL-97 and TW99 (just a marginal 
cover) is not satisfied either. The detailed criterions of the M-R constraints are presented in Table 
IVTl It is shown that the predictions given by DDRHF with PKO series and RMF with PKl, TMl, 



53l |. For comparison. 



DD-MEl, DD 
In Refs. 



ME2 and PKDD fulfill all the M-R constraints. 
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59(], the radius of neutron stars with the mass 1.4 Mq was found to be correlated 
with the neutron skin thickness of 208p]-, ^^yI as the symmetry energy. If the oberservation can 
limit the radius of neutron stars to a narrow range, a strong constraint can be imposed on the 
symmetry energy. On the other hand, if the neutron skin thickness of ^'^^Pb or the symmetry energy 
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FIG. 10: (Color online) Mass-radius relations of neutron stars provided by the DDRHF and RMF calculations 
and the corresponding maximum masses are marked in filled star symbols. For comparison, are also shown 
the four separate observational constraints from RX J1856 (gray grided region), 4U 0614+09 (cyan grided 
area), 4U 1636-536 (yellow grided area) and EXO 0748-676 (wine line for la error). The black region is 
excluded by causality that R > 2.9GM/c^ flQ- See the text for details. 



TABLE VI: The criterion of the M-R constraints: (1) the isolated neutron star RX J1856, (2) EXO 0748- 
676, (3) the low-mass X-ray binary 4U 0614-1-09, (4-u) 4U 1636-536 with its upper mass limits, and (4-1) 4U 
1636-536 with its lower mass limits. Fulfillment (violation) of a constraint is indicated with +(— ) and the 
marginal cover is marked with A. See the text for details. 
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could be precisely determined from the terrestrial experiments, it will be helpful to understand the 
neutron star structure and rule out some EoSs for the neutron star matter. As seen in Fig. [10] 
and Table IVl although several equations of state provide similar maximum masses, they still show 
some discrepancy for the radius of neutron stars with the mass 1.4 Mq. In Table IV] the DDRHF 
interactions predict this radius in a range from 13.8 km to 14.1 km, while the nonlinear RMF 
one NLSH gives the largest value 14.9 km, and the density dependent RMF one TW99 shows the 
smallest 12.4 km. All the calculated results except TW99 are coincident with the X-ray spectral 
analysis of the quiescent LMXB X7 in the globular cluster 47 Tuc, which requires a rather large 
radius of 14.5j[}'g km for 1.4 Mq compact stars {l?! . 
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IV. SUMMARY 

In this paper, the equations of state for symmetric nuclear matter, pure neutron matter and 
/3-stable neutron star matter have been studied within the density dependent relativistic Hartree- 
Fock (DDRHF) theory with PKOl, PK02 and PK03. Substantial effects from the Fock terms 
are found in describing the asymmetric nuclear matter at high densities. Due to the contributions 
from the Fock terms of a- and w-couplings, stronger density dependence on the symmetry energy 
is obtained from DDRHF at high densities, as compared to the RMF calculations with the density 
dependent meson-nucleon couplings. Because of the weak density dependence of gp in PKOl, 
which induces remarkable contributions from the p-meson coupling to the symmetry energy, PKOl 
shows stronger density dependence on the symmetry energy than both PK02 and PK03. With 
the obtained equations of state for /3-stable nuclear matter, the properties of neutron stars are 
investigated within the DDRHF theory for the first time and the recent observational constraints 
of compact stars are also introduced to test the applicability of the DDRHF models. 

Due to the extra enhancement from the a and lo exchange terms on the symmetry energy, 
large proton fractions in neutron stars are predicted by the DDRHF calculations, which affects 
essentially the cooling process of the star. For the DU process occurring, DDRHF with PK02 and 
PK03 gives critical neutron star mass ~ 1.45 Mq, which are close to the limit 1.5 Mq from the 
modern soft X-ray data analysis in the temperature-age diagram and fulfil the weaker constraint 
1.35 Mq. In contrast, fairly small mass limits are presented by the calculations of DDRHF with 
PKOl, RMF with the non- linear self-couplings of mesons, and RMF with PKDD, mainly due to 
their stronger p-coupling contributions to the symmetry energy at high densities. Different from 
these two cases, the RMF calculations with TW99, DD-MEl and DD-ME2 do not support the 
occurrence of DU process in neutron stars at all. In addition, the radii of 1.4 Mq neutron stars 
are correlated with the symmetry energy as well. In general, stronger density dependence on the 
symmetry energy leads to larger radius for 1.4 Mq neutron star. The radii given by the DDRHF 
calculations lie between 13.8 and 14.1 km, larger than the RMF calculations with the density- 
dependent meson-nucleon couplings, and smaller than the ones with the non-linear self-couplings 
of mesons except GL-97. 

For the maximum masses and central densities of neutron stars, they are tightly correlated with 
the behavior of the pressure with respect to the density. Due to the similar density dependent 
behaviors of the pressure, identical maximum masses (~ 2.5 Mq) of neutron stars are found in the 
calculations of DDRHF, and RMF with DD-MEl and DD-ME2, as weh as the central densities 
around 0.80 fm^^. The results are in reasonable agreement with high pulsar mass 2.08 it 0.19 Mq 
from PSR B1516-I-02B recently reported. The mass-radius relations of neutron stars determined 
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by the DDRHF calculations are also consistent with the observational data from thermal radiation 
measurement in the isolated neutron star RX J1856, QPOs frequency limits in LMXBs 4U 0614+09 
and 4U 1636-536, and the redshift limit determined in LMXBs EXO 0748-676, which are only 
partially satisfied in the RMF calculations with GL-97, NLl, NL3, NLSH and TW99. 
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